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For simple magnets quenched to below their critical 
point, one generally expects, in the scaling regime where 
t, s imicro an d also t— s 3> t m i CTO (i m icro is a microscopic 
reference time), the following scaling behaviour of the 
linear autoresponse function 
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where <f>(t) is the order parameter at the observation time 
t and h(s) is the magnetic field at the waiting time s. 

For some time, Corberi, Lippiello and Zannetti (CLZ) 
P, have advocated a phenomenological formula 
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where du,L are the upper and lower critical dimensions, 
respectively and n is the number of components of the 
order parameter. On the other hand, the commonly 
accepted physical picture of the ageing (coarsening) 
process going on after the quench has it that ordered 
domains form very rapidly and that the ageing comes 
from the movement of the domain walls between the or- 
dered regimes. If this idea is combined with scaling ar- 
guments and if one takes into account that in certain 
systems (such as the Ising model and referred to as class 
S) the spatial correlations decay exponentially while in 
others (such as the spherical model and referred to as 
class L) the spatial correlations decay algebraically, viz. 
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, one obtains 
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where z is the dynamical exponent. The result J3J) has 
been reproduced in many studies. However, if eq. @ 
of CLZ should turn out to be correct that would also 
invalidate the simple physical picture of the coarsening 
process mentioned above. A good practical way to decide 
between and the prediction © appears to be a study 
of the ID Ising model with a non-conserved order pa- 
rameter and quenched to T < T c from a fully disordered 
initial state. Here z — 2 is known [(| and from eqs. @ 
and |J3J one has a = 1/4 and a = 1/2, respectively. 



CLZ first arrived at eq. © by analyzing numerical 
data of the field-cooled susceptibility xzfc(^s) = 
f*du R(t,u). From a straightforward integration of QJ, 
they obtained xzFc(t, s) ~ s~ a and proceeded to extract 
a. However, it was pointed out that xzFc(ijs) does not 
obey a simple scaling but rather there may be a further 
and dominant contribution coming from the upper limit 
of integration. For systems of class S, there are well- 
defined domains with domain walls whose thickness is 
small with respect to the domain size and one has |5| 
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where xo,i are constants and Jm{x) > is a scaling func- 
tion. Indeed, it can be shown that A — z~ x — n 5], where 
the exponent k > describes the time-dependent scaling 
of the width of the domain walls w(t) ~ t K pfl. In the 2D 
Ising model k = 1/4 is known 0, hence A = 1/4 which 
explains the early result of CLZ but it becomes clear that 
the term naively expected from the scaling behaviour ifTjl 
merely yields a finite-time correction. CLZ did not take 
into account the condition i-s » imicro for the valid- 
ity of the scaling form (JIJ in their analysis and hence 
have missed the leading time-dependent term in ijlfl. 

On the other hand, for systems of class L, we argued 
previously that because of the long-range correlations, 
the effective width of the domains should be the same 
as their linear size which leads to A = |5(. Since the 
clusters have no 'inside', the term analogous to xo = in 
eq. JIJ is absent, but it's role is taken over by the constant 
Xi- CLZ [lj point out correctly that for both classes 
S and L the constant lim^oo XzFc(t, s) (with x = t/s 
fixed) and which equals xo for class S and Xo + Xi for 
class L can be related to the equilibrium magnetization 
through the static fluctuation-dissipation theorem. In 
summary one has, with a given by Q 
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In their comment, CLZ 0] now argue that the expo- 
nents A and a were really not distinct and that one rather 
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FIG. 1: (a) Scaling plot of the reduced field-cooled susceptibil- 
ity XzFc(i, s) — xo against t/s in the 2D Ising model quenched 
to T = 1.5, for waiting times s = 25, 50, 100, 200, 800, 1600 
(from bottom to top), (b) Comparison of the data with a fit 
-1 / 4 + ais" 1 ' 2 (black lines) and a fit 
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Xo = aos~ 1//4 (grey lines), for several values of 
t/s. For clarity the data for t/s = 2 and for t/s = 3 have 
been shifted upwards by a factor 1.4 and 1.2, respectively. 



had a — A. From now on, we concentrate on systems of 
class S and furthermore on the 2D Ising model quenched 
below T c . From numerical simulations, CLZ reproduce 



once more that xzfc(^, s) — xo 
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as generally 



expected. But they do not consider explicitly the correc- 
tions to that leading scaling behaviour and merely claim, 
but do not prove, the absence of the scaling corrections 
of order s" a , which are expected from eq. J5J. 

Indeed, the scaling behaviour of xzfc is not as simple 
as claimed by CLZ. We show this in figure QJi, with data 
obtained from the standard heat-bath method. It can 
be seen that the collapse in this scaling plot is far from 
complete which strongly indicates the presence of sizeable 
finite-time corrections. The origin of these becomes clear 
in figure^) where we compare our data to the asymptotic 
form xzfc(^7 s ) — Xo = »o s_1 ^ 4 + ais -1 / 2 , for several 
values of t/s. While taking into account both terms does 
reproduce the data well, we included for comparison also 
the fits where a± = was fixed by hand, as suggested by 
CLZ. It is clear that the data, which display a pronounced 
curvature, can only be fitted if ai < 0, as expected from 
©. Indeed, the best fits yield ao = 0.68 and a\ = —0.35 
for t/s = 2, a = 0.62 and a\ = -0.19 for t/s = 3 
and ao = 0.54 and a\ — —0.06 for t/s = 5 in excellent 
quantitative agreement with J5j. We conclude that there 
is no indication in favour of an absence of the last term in 
eq. (JHJ, hence a > A which provides clear counterevidence 
against the proposal of CLZ. 

We had already given earlier an analysis of the scaling 
of the thermoremanent magnetization (where terms of 
order s~ A do not occur). By taking into account also the 



leading finite-time correction we found full compatibility 
with a = 1/2 but inconsistency with a = 1/4 Q. 

We may also show that indeed a — 1/2 in the 2D 
Ising model through the so-called intermediate integrated 
response 
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which has the advantage that the leading term ~ t~ A 
which dominates over the scaling term in xzfc is absent. 
In figure |2 we examine the scaling of Xint(*> s ) an d try 
to achieve a collapse for a = 1/2 and a — 1/4. While 
there is a nice collapse for a = 1 /2 already for relatively 
small values of the waiting time (unless t/s is too close to 
unity), the data fail to collapse for a = 1/4. Again, this 
is fully consistent with a — 1/2 but excludes a = 1/4. 
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FIG. 2: Scaling plot of the intermediate susceptibility 
Xint it, s) for (a) the assumed value a = 1/2 (b) the assumed 
value a — 1 /4. Error bars are smaller than the symbol sizes. 

Turning to systems of class L, CLZ nicely reconfirmed 
the expected scaling forms eqs. ()■' '>!■")[) . However, there is 
nothing in their test with contradicts the simple scaling 
we used in [5j to obtain A = 0. 

In conclusion, having reexamined the scaling of some 
integrated susceptibilites, we have shown that the scaling 
of Xzfc (A s ) does indeed contain at least two important 
contributions, see eq. JSJ, which is against the proposal 
of CLZ. We stress that a simple demonstration of scaling 
of Xzfc is not enough to be able to reliably know which 
of the exponents A or a is measured. It is more safe to 
study a quantity such as the intermediate susceptibility, 
which does not suffer from this difficulty. Applied to the 
Ising model quenched to T < T c , our results fully confirm 
eq. (El but disagree with eq. (J2J- 
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